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In this paper, a Hirota method is developed for applying to the nonlinear 
Schrodinger equation with arbitrary time-dependent linear potential which denotes 
the dynamics of soliton solutions in quasi-one-dimensional Bose-Einstein condensa- 
tion. The nonlinear Schrodinger equation is decoupled to two equations carefully. 
With a reasonable assumption the one- and two-soliton solutions are constructed 
analytically in the presence of an arbitrary time-dependent linear potential. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) jl, 2| which strongly stimulates the 
exploration of nonlinear properties of matter waves have opened a new field of nonlinear 
atom optics, such as four wave mixing in BEC's the study of various types of excita- 
tions. One of particular interest is macroscopically excited Bose-Einstein condensed states. 



such as vortices and solitons |7H13|. The existence of solitonic solutions is a general 
feature of nonlinear wave equations. For the case of an atomic Bose-Einstein condensate, 
the macroscopic wave function of the condensate obeys the so-called Gross-Pitaevskii (G- 
P) equation, whose nonlinearity result from the interatomic interactions. It is well known 
that the G-P equation has of either dark or bright solitons depending on the repulsive or 



attractive nature of the interatomic interactions, respectively. A dark soliton |l4l Il5| in 
BEC is a macroscopic excitation of the condensate which is characterized by a local density 
minimum and a phase gradient of the wave function at the position of the minimum. A 
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bright soliton 16Nl8| in BEC is expected for the balance between the dispersion and the 
attractive mean-field energy. Several methods have been applied to obtain the soliton so- 
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20j, as well as the dynamics of 



lutions of G-P equation with different potential |9l-ll3. 
the excitation of the condensate was discussed. When the longitudinal dimension of the 
BEC is much longer than its transverse dimensions which is the order of its healing length, 
the G-P equation can be reduced to the quasi-one-dimensional (quasi-lD) regime. This 



trapped quasi-low-dimensional 



21| condensates has offered an useful tool for investigating 



the nonlinear excitations such as solitons and vortices, which are more stable than in 3D, 
where the solitons suffer from the transverse instability and the vortices can bend. Thus 
the study of both theory and experiment is very important for the soliton excitations in 
quasi-low-dimensional BECs. 

In this paper, we consider the mean-field model of a quasi- ID BEC trapped in a linear 
time-dependent potential which is given by 

where J dX = N is the number of atoms in the condensate, /(T) is the arbitrary 
function of time T, the interacting constant of two-atom is given hj g = 2h?a/mll ( 22| |) 
with m the mass of the atom, a the s-wave scattering length (a > for repulsive interaction; 
while a < for attractive interaction), and /q = ^yh/muo denoting the characteristic length 
extension of the ground state wave function of harmonic oscillator. Making a dimensionless 
transformation, we can rewrite equation (1) as 

+ + x/ (t) 7^ + IV^I V = 0, (2) 

where x is measured in units of /q, ^ in units of mt^jh^ ifj in units of the square root of A^/q, 
the interaction constant n is defined as yU = —2NlQa, and / (t) = ~-j^f{j-j^)- The exact 
soliton solutions of Eg. (121) can been constructed by the inverse scattering method [23| and 
F-expansion method [2^]. It should be noted that in the absence of the linear potential, i.e. 
x/ (t) = 0, a Hirota method can be applied to Eq. directly for getting the bright and 
dark soliton solutions. However, with the consideration of the linear potential the Hirota 
method should be developed carefully. This is our purpose in the present paper. With 
a reasonable assumption we demonstrate how to construct the exact one- and two-soliton 
solutions of Eq. (12]) in terms of this developed technique. 
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II. ONE-SOLITON SOLUTION 

Now we introduce the main idea of the Hirota method briefly. Firstly, it apply a direct 
transformation to the nonlinear equation. Then in terms of the reasonable assumption 
the nonlinear equation can be decoupled to two equations from which the one- and two- 
soliton solutions can be constructed effectively. To this purpose we consider the following 
transformation 

where G {x, t) is complex function and F (x, t) is a real function. With this transformation 
Eq. (|2]) becomes 

F (^Dt + ]^Dl^ G-F + G-F^- xf{t) - G{^DIF ■ F - iiGG) = 0, (4) 

where the overbar denotes the complex conjugate, Dt and are called the Hirota bilinear 
operators defined as 

D-D-G ix, t) . F {X, t) = (^^ - ^ J - ^ j G (X, t) F (x', t') |_,,^, . (5) 

If the term xf{t) = Eq. (j4]) reduces to the normal nonlinear equation which can be 
decoupled easily to two equations. In the presence of the term xf(t) we should deal with 
Eq. (jl]) carefully. Many attempts show that Equation (jl]) can be decoupled as 

^DlF .F-f,GG = 0, 
iDt + ^Dl + xfit)^G-F = 0, (6) 

in which the spatial and time dependence term xf{t) will give a difficulty for getting solutions 
as shown below. Now the Eq. (jS]) has made the Eq. ([2]) to the normal procedure of Hirota 
method for getting the exact soliton solutions. By making a series of suitable assumption 
for the expression of G and F, the exact one- and two-soliton solution can be obtained 
analytically. In order to obtain the bright one-soliton solution of Eq. ([2]) which correspond 
to the case a < 0, we proceed in the standard assumption 

G = xGu F=l + x'Fu (7) 
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where x is an arbitrary parameter which will be absorbed in expressing the soliton solution 
in the following sections. Substituting Eq. ([7]) into Eq. ([H]), then collecting the coefficients 
with same power in x, we have 
(1) for the coefficient of x 

^Dt + \Dl + xf{t)]G^■l = Q, (8) 



2 

(2) for the coefficient of x^ 

DlFi - fiGiGi = 0, (9) 

(3) for the coefficient of x^ 

^D, + ^Dl + xf{t)^G^■F^=0, (10) 

(4) for the coefficient of 

DlF, ■ Fi = 0. (11) 

Using the definition ([5]) the above equations can be expressed in detail. For example, in 
order to satisfy Eq. we can assume Gi has the form 

Gi = e''i. (12) 
Substituting Eq. (IT^ into Eq. we have 

iVit + \]\xx + \^\x + ^/(^) = 0. (13) 

Because of the presence of the term xjif) there are some difficulties for getting general 
solutions of the above equation. For some conveniences in this paper we assume t]^ has the 
form 

r/i = Pi(t)x + f]i(t), (14) 

with the time- dependent functions -Pi(t) and fii(t) to be determined. With the restriction, 
i.e., Eqs. ^ and ( ITT]) we get 

Fi = exp(r7i+r7i + Aii), (15) 

where 

/ln = ln ^ 



(Pi + Pi) 



2 • 
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Substituting the solutions (fT^ and (fT5l) into Eq. (fTU]) we obtain the equations of Pi(t) and 
fii(t) as 

[^Pi,i(t) + fit)] X + + ^P^ = 0, (16) 

form which we can determine the expression of Pi{t) and Qi{t). It is obvious that Eq. f lT6|l 
imphes the natural conditions 



^Pi,i(t) + fit) = 0, 
zfii^t) + ^P'(t) = 0. 
Solving the above two equations, we get the expression of Pi and Qi as 

Pi = z f f{r)dr + U, 
Jo 

j2 



(17) 



(r) rfr + C 
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(18) 







where C,io and (iq are complex parameters in general. With the Eqs. f lT2|) . f lTSj) and 
after absorbing x, the bright one-soliton solution of Eq. (|2]) can be derived as 



I _|_ g»?i+-'yi+Aii ' 



(19) 



where 



Au = In 



t 

I I /(r)rfr + ^10 





+ lfpHr)dT + C,,. 



T \ 2 • 



(20) 



(^lO + U)' 

In the case of f{t) = 0, one-soliton solution f[T^ can reduce to the solutions of the normal 
nonlinear Schrodinger equation. When f{t) = constant the solution (1191) is the same results 
( 2^) reported earlier. As f{t) = bi + lcos{iot), the solution fl5^ is the same results ( 24|). 
From Eq. fl20l) we can see the linear time-dependent potential can change the soliton velocity 
and frequency. 



III. TWO-SOLITON SOLUTION 



In this section we will give the analytical expression of two soliton solution of Eq. (I2j). 
To this purpose we now assume that 

G = xGi + x'G2, F = l + x'Fi + x'Ps. (21) 
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By employing the same procedure before we obtain the following set of equations from Eq. 
dn]), corresponding to the different powers of x 

(1) for the coefficient of x 

^A + ^I^^ + x/(t))Gi-l = 0, (22) 

(2) for the coefficient of x^ 

DlFi ■ 1 = fiGiGi, (23) 

(3) for the coefficient of x^ 
(^iDt + ^Dl^ (Gi ■ Fi + • 1) + xf{t) {GiF, + G2) = 0, (24) 

(4) for the coefficient of 
DlFi ■ Fi + 2DIF2 ■ 1 = 2fi{GiG2 + G2G1), (25) 

(5) for the coefficient of x^ 
(^iDt + ^Dl^ (Gi ■F2 + G2- Fi) + xf{t) (G1F2 + G2F1) = 0, (26) 

(6) for the coefficient of x^ 

DlFi ■ F2 = fiG2G2, (27) 

(7) for the coefficient of x'' 

+ ^Dl + xf{t)^ G2-F2 = 0, (28) 

(8) for the coefficient of x^ 

DIF2 ■ F2 = 0, (29) 

As discussed in the one-soliton solution, we can solve the equations from (1221) to (12^ in turn 
for getting the expression of G and F. In order to construct the two-soliton solution of ([2]) 
we assume Gi has the form 

Gi = exp Tji + exp 7^2, (30) 

where 

Vj = Pj{t)x + Qj{t), 
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in which the time- dependent functions Pj(t) and flj(t), j = 1,2, to be determined. Substi- 
tuting Eq. fl^ into the relation (1221) we have 

I [iPi^t + fit)] X + iQi^t + ^P^l + I [zPs.t + f{t)] x + ifl^.t + | exp r^^ = 0, 



which implies that 

^PJ,t + fit) = 0, 

+ = 0. 

where j = 1,2. From the above equations one can find the solutions 

Pj it) = ^ f f{r)dr + 
Jo 

where and (j^, j = 1,2, are complex parameters in general. Combining Eq. fl30|) with 
Eq. (l23i) we obtain the expression of Fi as 

Fi = exp(r7i + r^i + An) + exp(r72 + + ^21) 

+ exp(r7i + r/2 + A12) + exp(r72 + r/2 + ^22), (31) 

where 

= In ^_ n , m, = 1, 2. 

+P 1 

With the help of Eqs. ([30]) and ([3ID we can simplify Eq. as 

iG2t + ^Ga^^ + xf{t)G2 

/i (P2 - Pl)^^^^^+^^+^^ ^ /i(^2-Pl)^_ ^r;i+fi,+r,,^ ^33) 



(Pi + Pi) (P2 + Pi) (Pi + P2) (P2 + P2) 

which shows that the expression of G2 has the form 

G2 = exp(?7i + V1 + V2 + ^1) + exp(r7i + ^72 + ^72 + ^2), (33) 
where the parameter 5j, j = 1, 2, is given by 

5, = lnf^yj = l,2, (34) 
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7i = (Pi + Pi) (P2 + Pi) , 72 = {Pi + ^2) (P2 + P2) 

/i(P2-Pl)' p /i(P2-Pl)' 



(Pi + Pi) (P2 + Pi) ' (Pi + P2) (P2 + P2) ■ 

Now we have obtained the expression of Gi, G2, and Pi in Eq. ([21]). Substituting Eqs. fl5U]) . 
Plj) and (ESD into Eq. we obtain 

F2 = e^i+^i+''2+%+«^ (35) 

where 

, /i^iP2-Pir 

K = In 



(Pl+Pl)^(P2 + P2) |Pl + P2r 

With the help of Eqs. ([30]), ([3l]), ([33]), and ([35]) one can find the Eqs. ([26]) and ([29]) are 
satisfied to the moment after a tedious calculation. So the two-soliton solutions has the form 

= J. (36) 

where 

G = e^^ + -|- g'7i+'7i+^2+'5i _|_ g'7i+%+f?2+'52^ 

_P — ]_ _|_ _j_ gr?2+%+-42i _j_ gr?i+rj2+Ai2 _|_ gr?2+rj2+A22 _|_ qVi+Vi+V2+V2+i^ 



Jo 

n,{t) = Pf (t) dr + c 



Sj = In 



P, 



where j = 1,2, and 



71 = (^10 + ^10) (^20 + "^lo) 

72 = (^10 + ?2o) (^20 + ^20) 

_ (^20 ~ Cio) 

(eio + ^io) (^20 + ^10) 

_ (^20 ~ Clo) 

i^w + ^20) ('C20 + ^20) 
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/^^ 1^20 ^10 1 



(■ClO + 'Clo) ('C20 + 'C20) |'ClO + ^2o| 



^mn = In ^_ 2 ' m, n = 1, 2. 

When f{t) = 0, the solution f l36|) denotes two sohton interaction of the normal nonlin- 
ear Schrodinger equation. The new expression (136|) implies that Hirota method has more 
advantage for getting new soliton solutions as well. 



IV. CONCLUSION 



In this paper, we investigate the soliton solutions of the nonlinear Schrodinger equation 
with an arbitrary time-dependent linear potential which denotes the dynamics of quasi-one- 
dimensional Bose-Einstein condensation. A developed Hirota method is applied carefully to 
the nonlinear Schrodinger equation. In terms of this developed technique we decoupled the 
nonlinear Schrodinger equation into two equations. Moreover, with a reasonable assumption 
the exact one- and new two-soliton solutions are constructed effectively. Our soliton inter- 
action will have useful application in the studies of Bose-Einstein condensation and optics 
communication in which the nonlinear Schrodinger equation are used widely. 
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